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Future hadron collider experiments will allow for a detailed investigation of ti spin correlation effects. In this 
talk, recent progress in the theoretical description of these effects is reported. In particular, next-to-leading order 
results for the ti spin correlations in qq — ► tiX are presented, using various spin quantization axes. 



Top quark pair production and decay can be 
described by perturbative methods. This is so 
because the top quark width is large compared to 
the QCD hadronisation scale and therefore the 
top decays before hadronisation takes place. The 
full information of the top quark interactions in- 
cluding spin correlations between t and i can thus 
in principle be obtained from studying differential 
distributions of the decay products. In this talk, 
we will report on recent progress in the theoretical 
description of top quark spin correlation effects at 
hadron colliders. 

Consider the reactions 

( 21 + n > 2 jets + P™ iss 
pp,pp -> tiX -> \ I + n > 4 jets + P£ liss (1) 
[ n > 6 jets. 

To leading order in the QCD coupling, two parton 
reactions contribute: 

qq, gg — + tt — > bW + bW~ — > 6 fermions. (2) 

The amplitudes for these two- to six-body pro- 
cesses, with intermediate top quarks of non-zero 
total width, were given first in [ [jj. The calcu- 
lation of the fully differential cross sections for 
the reactions ([!]) including NLO QCD correc- 
tions is a formidable task. A substantial simplifi- 
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cation in the computation of the radiative cor- 
rections is achieved if one performs an expan- 
sion in Tt /mt and Tw / fnw and keeps only the 
terms with the highest degree of resonance. In 
this leading pole approximation [ ^ the radiative 
corrections can be classified into so-called fac- 
torisable and non-factorisable corrections. The 
non-factorisable NLO QCD corrections were cal- 
culated in [ ID . For the factorisable corrections 
the square of the complete matrix element A4^' 
for the partonic reactions (|^) is of the form 

| M^\ 2 (x Tr [pR^p] = Pa > a R { 3,f30'P0'0- (3) 

Here denotes the density matrix for the pro- 
duction of on-shell tt pairs, the label A indicates 
the process, and /?, p are the density matrices de- 
scribing the decay of polarised t and t quarks, 
respectively, into specific final states. The sub- 
scripts in (H) denote the t, i spin indices. Note 
that both the production and decay density ma- 
trices are gauge invariant. 

For the process qq — > ti, which is the dominant 
process at the Tevatron, the production density 
matrix R is given in terms of transition matrix 
elements as follows: 

p (t a tp\T\qq)(qq\Ti\t a >ip>) 
Haa',/30' = M ' ^ 

colors JV 99 
qq spins 

where the factor N q q = (2N C ) 2 = 36 averages 
over the spins and colors of the initial qq pair. 
The matrix structure of R is 



Raa',/3/3' = A5 aa > 5pp> + C t j (V ) aa ' ((7° 

+ Bi((j l ) aa ,8f3p, + Bi5 aal (a l )pp> ,(5) 

where a % are the Pauli matrices. Using rota- 
tional invariance the 'structure functions' Bi,Bi 
and Cij can be further decomposed. The function 
A, which determines the it cross section, is known 
to next-to-leading order in a s from the work of [ 
0, ||. The corresponding function for the gluon 
fusion process is also known in NLO [ |], |(|. Be- 
cause of parity (P) invariance the vectors B,B 
can have, within QCD, only a component normal 
to the scattering plane. This component, which 
amounts to a normal polarisation of the t and t 
quarks, is induced by the absorptive part of the 
scattering amplitude, and it was computed for qq 
and gg initial states in [ |?], || to order The 
normal polarisation is quite small, both for it pro- 
duction at the Tevatron and at the LHC. Parity 
and CP invariance of QCD dictates that the func- 
tions C%j, which encode the correlation between 
the t and i spins, have the structure [ |^] 

Qj = ciSij + c 2 piPj + c 3 kikj + C4,(kifij +pikj), (6) 

where p and k are the directions of flight of the 
initial quark and of the t quark, respectively, in 
the parton cm. frame. The production density 
matrix for the reaction gg — > ti can be defined 
and decomposed in an analogous fashion. To 
Born approximation the functions c r were given, 
e.g., in [ Theoretical studies of spin correla- 
tions have been performed at leading order in a s 
in [ [Icj, [yj H II II HI- An attempt to detect 
spin correlations in a small ti dilepton sample col- 
lected at the Tevatron was recently reported by 
the DO collaboration [ |l6) . A feasibility study for 
the LHC can be found in [ |l7||. 

We have recently computed the ti spin density 
matrices for the parton reactions qq — > tt, tig to 
order [|lS[]. We have also calculated, for these 
reactions, the degree of the ti spin correlation at 
NLO for different t and i spin quantization axes. 

The ti spin correlations can be inferred from 
appropriate angular correlations and distribu- 
tions of the t and i decay products. In the SM the 
main top decay modes are t — > bW — > bqq' , blvi. 



Among these final states the charged leptons, or 
the jets from quarks of weak isospin -1/2 originat- 
ing from W decay, are most sensitive to the po- 
larisation of the top quarks. The one-loop QCD 
corrections to the semileptonic decays of polarised 
top quarks and to t — > W + b can be extracted 
from the results of [ [lTJ] and [ |2(| , respectively. 
In the following we describe our computation of 
the density matrices for ti production by qq an- 
nihilation. At NLO we have to consider the reac- 
tions 



q{Pi) + q(P2) 
and 



t(h)+t(k 2 ), 



q(pi) + q(p2) -> Kh) + i(k 2 ) + g (k 3 ) 



(7) 



(8) 



In order to determine these functions to order 
a^. we first computed the one-loop diagrams that 
contribute to (|j) . Dimensional regularization was 
employed to treat both the ultraviolet and the in- 
frared and collinear singularities which appear in 
the diagrams. The ultraviolet singularities were 
removed by using the MS prescription for the 
QCD coupling a s and the on-shell definition of 
the top mass mt- The initial quarks are taken 
to be massless. After renormalisation the density 
matrix for the ti final state still contains single 
and double poles in e = (4 — D)/2 due to soft 
and collinear singularities. These poles are can- 
celled after including the contributions of the re- 
action (H) and mass factorization. For the latter 
we used the MS factorization scheme. We avoided 
the computation of the exact density matrix for 
the reaction (^|) in D dimensions by employing a 
simple version of the phase-space slicing method 
We divided the phase space into four re- 
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gions, namely the region where the gluon is soft, 
the two regions where the gluon is collinear to 
one of the initial state massless quarks (but not 
soft), and the complement of these three regions, 
where all partons are 'resolved'. This decomposi- 
tion can be performed using a single dimension- 
less cut parameter ir m i n . For example, the soft re- 
gion is defined by the requirement that the scaled 
gluon energy in the cm. system x g = 2E g /^/s is 



smaller than x min . In the soft region we used 
the eikonal approximation of the matrix element 
for reaction (H) and the soft limit of the phase 



space measure. The integration over the gluon 
momentum can then be carried out analytically 
in D dimensions. The two collinear regions are 
defined by {cos9 qg > (1 — x m - ln ) and x g > x m - ul ) 
and (cos8 qg < (— 1 + x min ) and x g > x min ), re- 
spectively, where 9 qg is the angle between the 
gluon and the quark in the qq cm. frame. In 
these regions we used the collinear approxima- 
tions for both the squared matrix element and 
the phase space in D dimensions. Finally, the ex- 
act spin density matrix for reaction (||) in four 
space-time dimensions was used in the resolved 
region, where all necessary phase space integra- 
tions can be carried out numerically. By construc- 
tion, all four individual contributions depend log- 
arithmically on the slicing parameter x m in, but 
in the sum only a residual linear dependence on 
£min remains, which is due to the approximations 
made in the soft and collinear regions. By varying 
£min between 10 -3 and 10 -8 we checked that for 
£min < 1CP 4 this residual dependence is smaller 
than our numerical error (which is less than a 
pcrmill for all results discussed below). 

After mass factorisation we are left with finite 
density matrices for the ti and the ti + hard gluon 
final states. As a check of our calculation we first 
compute the total cross section for qq — > ti + X 
at NLO. If one identifies the MS renormalisation 
scale fi with the mass factorisation scale [If and 
neglects all quark masses except for m t , then one 
can express the cross section in terms of dimen- 
sionless scaling functions [ H : 



r(s,m?) 



(9) 



where s is the parton cm. energy squared and 
i] = s/Arni — 1. We have compared our result for 
ct qq as a function of rj with those of [ [|, |5| and 
found perfect agreement. 

We now consider the following set of spin- 
correlation observables: 



o 4 = 4(p*. Sl )(pr -s 2 ), 

06=4(di-si)(d2-s 2 ), 



(13) 
(14) 



where Si,S2 are the t and i spin operators, re- 
spectively. The factor of 4 is conventional. With 
this normalization, the expectation value of 0\ is 
equal to 1 at the Born level. The expectation val- 
ues of the observables (0), (§J), ©, and (0) 
determine the correlation of different t, i spin pro- 
jections. Eq. ( |i"T|) corresponds to a correlation of 
the t and t spins in the helicity basis, while jl^ ) 
correlates the spins projected along the beam line 
in the parton c.m.s. The 'beam-line basis' used 
T3~|) was defined in [ 13 and refers to spin axes 



being parallel to the antiquark direction in the t 
rest frame p 2 and to the quark direction in the i 
rest frame f>\* 1 respectively. The spin axes di 2 
i n (|ll|) c orrespond to the so-called 'optimal basis' 
[ [23|, |14[ to be discussed below. 

For quark-antiquark annihilation it turns out 
that the spin correlation (|l2|) [ [To| , |l5[ and the 
correlation in the beam-line basis ([l 3|) [ |l3| are 
stronger than the correlation in the helicity ba- 
sis. A spin-quantization axis was constructed in 
[ |23|, Q with respect to which the t and i spins 
are 100% correlated to leading order in the QCD 
coupling, for all energies and scattering angles. In 
terms of the structure functions of (jg) this means 
that the 'optimal' spin axis d fulfills the condition 



diCijdj — A. 



(15) 



The existence of a solution of Eq. (|15|) is a special 
property of the leading order spin density matrix 
for the reaction qq — > tt. One finds [ 0, |l4| : 



-pi + (1 - 7i)(Pi ' ki)ki 



l-(px- k!)a(l-7?) 



(16) 



Oi = 4si • s 2 , 



(10) 



where 71 = Ei/m t . The construction of this axis 
explicitly uses the leading order result for the spin 
density matrix, and different generalizations to 
higher orders are possible. We use in ( |l4| ) as spin 
axes: 



2 =4(k 1 -s 1 )(k 2 -s 2 ), 
3 = 4(pi • Si)(pi • s 2 ), 



(11) 
(12) 



-pi + (1 - 72) (pi • k 2 )k 2 



i-( P1 .k 2 ) 2 (i-7i) 



(17) 



where 72 = B^/rnt- For the 2 to 2 process qq — ► 
tt, d 2 = di = d. 

The expectation value of a spin-correlation ob- 
servable O at parton level can be written at next- 
to-leading order in analogy to (||) as follows: 

{0) q ~ q = <$(»»)+ 47^0$ fa) 

+ S#fa)M/4M a )). (18) 

Note that these quantities depend explicitly on 
the factorization scale fip, but only implicitly 
(through a s ) on the renormalization scale \x. This 
is because a factor a 2 s drops out in the expecta- 
tion values, and hence the Born result is of order 

Our results for the functions gW{r)), g^iv) 
and ffp^ are shown for the five observables 
(jjj)~dll) m Figs. 1-5. In each figure, the dot- 
ted line is the Born result gg^(rf), the full line 
shows the function g q q (77) , and the dashed line is 

g qq - (rj). A general feature of all results is that 
the QCD corrections are very small for values 
of t) < 1. For larger values of 77, the functions 
3^(77) depart significantly from zero. Also, the 

functions g qq (77) become nonzero, with less dra- 
matic growth as 77 — > 00 and with an opposite 
sign as compared to g q ^ (77) . The phenomenologi- 
cal implications of these features for spin correla- 
tions at the Tevatron and the LHC will be studied 
in detail in a future work. Here we merely note 
that the substantial QCD corrections for large 77 
will be damped by the parton distribution func- 
tions, which decrease rapidly with 77. Moreover, 
at Tevatron energies values of 77 above ~ 30 are 
kinematically excluded. 

To summarize: We have computed the spin 
density matrices describing tt production by qq 
annihilation to order q^. Further we have evalu- 
ated the scaling functions encoding the QCD cor- 
rections to spin correlations, using a number of 
different spin quantization axes. This work pro- 
vides a building block, which was missing so far, 
towards a complete description of the hadronic 
production of top quark pairs at NLO in the 
strong coupling. 
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Figure 1. Dimensionless scaling functions g qq ^(r)) 

(dotted), g^(ri) (full), and £fg^ (77) (dashed) that 
determine the expectation value {0\) qq . 
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Fi gure 2. Same as Fig.l, but for (fD2) qq . 
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Figure 3. Same as Fig.l, but for (0 3 } qq . 
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Figure 4. Same as Fig.l, but for {Ot±) q q. 




Figure 5. Same as Fig.l, but for (O5 
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Discussion 

A. P. Contogouris, Univ. of Athens 
I am worried not by the method of your calcu- 
lation, but by the importance of your numerical 
results. In particular the channel gg — > ti(g) can 
change them very much, and I am glad that you 
have already set up its calculation. Of course, 
qq — > tt(g) is required to be complete. 
A. Brandenburg 

/ agree. Note that the plots show results at the 
parton level for qq — > tt(g). 



